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We study deeply virtual Compton scattering off a spin-one particle, as it is the case for the coherent 
scattering off a deuteron target. We extend our approach, formulated initially for a spinless case, 
and discuss the role of twist three contributions for restoring the gauge invariance of the amplitude. 
Using twist three contributions and relations which emanate from the QCD equations of motion, 
we derive the gauge invariant amplitude for deeply virtual Compton scattering (DVCS) off hadrons 
with spin 1. Using the derived gauge invariant amplitude, the single spin asymmetry is discussed. 

PACS numbers: 12.38.Bx, 13.60.Le 



Deeply virtual Compton scattering (DVCS) off the deuteron target has recently attracted much attention 
from the experimental point of view UM- One of the main reasons of this interest is the fact that the DVCS 
process gives information about the generalized parton distributions (GPDs). From the theoretical point of 
view, the leading twist-2 GPDs for the deuteron were studied in [5|-u|. However, the leading twist- 2 accuracy 
for the DVCS amplitude, calculated in thecase where the final deuteron gets a significant transverse momentum, 
is not enough for the study of such processes. This is mainly due to the fact that there is a problem with the 
photon gauge invariance of the DVCS amplitude in the leading twist-2 approximation in the Bjorken limit (see, 
for example, |8|). Besides, the relevant terms are proportional to the transverse component of the momentum 
transfer and provide the leading contribution to some observables. 

Extending the Ellis-Furmanski-Petronzio-Efremov-Teryaev (EFPET) approach (see, 0, [l(|) to the non- 
forward case, this problem was first solved in [Tl| where it was demonstrated how the inclusion of twist-3 
contributions related to the matrix elements of quark-antiquark-gluon operators, can restore the gau ge i nvari- 
ance of the DVCS amplitude off a (pseudo) scalar particle (e.g. pion, He 4 ). Then, the main ideas of [ljj were 
used and generalized for the nucleon target (see, jl2Hl7j and [18[ for different processes). 

In this paper, we adhere to the approach [TlJ and make a comprehensive analysis of the twist three contribu- 
tions to the DVCS amplitude off a spin-1 hadron 1 . Since the parameterizations of the relevant hadronic matrix 
elements of th e q uark-gluon operators depend much on the spin of the external hadrons, following the very 
useful idea of [121 ] we start our study with a parametrization-free approach, and then apply it to the specific 
case of spin-1 hadron. 



I. V. Anikin/'^B R. S. Pasechnik, 3 '0 B. Pire, 4 <0 and O. V. Teryaev 2 '! 



I. INTRODUCTION 




Actually, the method described in and in this paper is suitable for a study of DVCS off an arbitrary spin hadron. 
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II. KINEMATICS AND APPROXIMATIONS 

Let us start with the discussion of the kinematics and approximations which we use in this paper. The process 
we consider is 

1 *(q)+D(p 1 )^ 1 (q') + D(p 2 ). (1) 

Here, we mainly focus on the deuteron as a target but all our approach is suitable for any spin-one hadron target. 
At the Born level, the Feynman diagrams corresponding to the considered process are depicted in Fig. [TJ This 
process is a hard exclusive reaction for which a QCD factorization theorem applies. In this case, the virtuality 
of the initial off-shell photon is usedas the large scale, i.e. q 2 — —Q 2 — > oo, while the final photon is on-shell 
with q' 2 = 0. Besides, this asymptotical regime is identical to the light-cone formalism. Therefore, we first 
introduce a light-cone basis which is constructed by the "plus" and "minus" vectors: 

n* = A(l, 0, 0, 1), n = -L(l, 0,0,-1), n* ■ n = 1 , (2) 

where A is an arbitrary and dimensionful constant which can be expressed via the Lorentz invariants. The exact 
form of A as a function of invariants depends on the frame which one works in. 

In the present paper, we consider the DVCS amplitude up to the twist three accuracy, discarding the con- 
tributions associated with the twist four and higher. Such a constraint imposes the following relations for the 
hadron relative and transfer momenta : 

P = = n H n ~ n , 

2 2 

A = p 2 - Pi = -2(P + 2(M 2 n + A T « -2£P + A T , 

P-A = 0, A 2 =t«0. (3) 

Notice that keeping the M 2 -term in the Sudakov decomposition of the relative momentum P (see, ([3])) leads to 
the necessity to include the twist four contributions as well, which goes beyond the scope of the present paper. 
Since corrections of the order 0(A|, / 'Q 2 ) demand a special care, at this moment, we postpone the study of them 
until a forthcoming paper. Notice that the detailed analysis of these contributions has recently been presented 
in P1I23. " " 

It is also instructive to introduce the photon relative momentum: 

Q = ^=q-J = q'+J, q' = (P-q')n, (P ■ q) = (P ■ Q) = (P ■ q) . (4) 

One has to emphasize that the approximations discussed in this section do not affect the generality of our study 
and can be applied to a study of arbitrary spin hadrons. 

III. FACTORIZATION AND THE GAUGE INVARIANT AMPLITUDE 

In this section, we briefly describe the factorization procedure applied to the DVCS amplitude up to the twist 
three accuracy. The details of this factorization can be found in [l8|, [2l|, [22[ . 

At the Born level, the sum of the amplitudes with the two- and three-particle correlators, or T^J and T$ 
amplitudes, has the following form (for the DIS case, see [Iol|): 



T$+T$ = J d 4 ktr E^(k)T{k) + J d 4 h d 4 k 2 tr 



E t _ Lpu {k 1 ,k 2 )T p (k 1 ,k 2 ) 



(5) 



where and E^ pu are the coefficient functions, at the Born approximation, with two quark legs and two 
quark and one gluon legs, respectively. In Eq. ((5]), we use the following notations: 



Tagik) = - J d A ze« k - A /V z (p 2 \^(z)M0)\pi), 

— OO 

r£fl(k,fe>) = - / d i z 1 d 4 z 2 e^- A / 2 ^+^~ k ^ ( P2 \^(z 1 )gA"(z 2 )M0)\Pi) > (6) 





Figure 1: The Feynman diagrams corresponding to the deeply virtual Compton scattering. The DVCS amplitude 
with two-particle correlators is depicted on the left panel, while the amplitude with the three-particle correlators - 
on the right panel. Notations: P = pi, P' = p 2 , K = k - A/2 w xP - A/2, K' = k + A/2 » xP + A/2, 
L = fci — A/2 ~ iiP — A/2, L = k2 + A/2 ~ a?2-P + A/2. Here, k and ki correspond to the loop momenta in the 
diagrams. 

where the underlined indices a, (3 denote the Dirac spin indices while the other indices correspond to the Lorentz 
ones. 

It is convenient to choose the axial gauge condition for gluons, i.e. n ■ A = 0, where n is the light-cone vector 
defined in Eq. Then, we carry out a decomposition of the loop momentum k over the light-cone vectors 
(the Sudakov decomposition) as follows 



k ll = xP p + {k ■ P)n M + fcj « x P M + fcj , x = k-n. 



(7) 



As the next step of the factorization procedure, we perform the replacement for the integration momentum in 
Eq. (0 as 



d A ki — > d A ki dxi5(xi — ki ■ n). 
This allows us to expand the two-quark coefficient function (see, (j^Jl) in a Taylor series: 

dE^(k) 



E^(k) = E^(xP) 



dk a 



(k-xP) a + ..., (k-xP) a ntki, 



k=xP 



Then, using the collinear Ward identity (see, fioj ) 

dkP 



(k,k). 



we, finally, arrive at the factorized (in the momentum space) DVCS amplitude which reads 



T$+T$= I dxtr 



where w pp < = S pp > — n p rP p , and 



E^(xP)T(x) 



+ J dx\ dx2 tr 
-l 



E flpv (xx,X2)Upp>T p i(xi,X2) 



(8) 



(9) 



(10) 



(11) 



To^(x) = - 

T p af3 (xi,x 2 ) 



dX 
2^ 



(p 2 |V'a(An)V' / 3(0)|pi), 



d\i d\2 
~2n ~2tT 



(pa|Va(Am) %DP (A 2 n)^(0)|pi), 



(12) 



where i_D M =i<9 M +gA^ is the QCD covariant derivative in the fundamental representation. This amplitude is 
also needed to be "factorized" in the Dirac space. This can be reached by making use of the Fierz decomposition 
over spinor indices. 

(2) 

In fact, the contributions of X)V are not completely independent from one another because of the QCD 
equations of motion (e.o.m.) for fermions. The next step is to use the QCD e.o.m. in order to reexpress 
the contributions of the correlators with the covariant derivative through the correlators which include only ip 

(2) 

and ip fields. Afterwards, the contribution of T^J presented in terms of the two-particle correlators should be 
combined together with the contribution of T$ in order to get the gauge invariant DVCS amplitude at the 
Born level. 

Let us now focus on the QCD equations of motion. For the sake of simplicity, we start within the approxima- 
tion where the three-particle correlators are absent. Below, we will discuss the extension up to the dynamical 
(or genuine) twist contributions related to the three-particle correlators. However, to derive the gauge invari- 
ant amplitude, it is sufficient to consider only the kinematical twist contributions since the kinematical and 
dynamical twists enter in the QCD e.o.m. and the DVCS amplitude additively (see [H EE El HI)- 

Let us consider the e.o.m. in the following form: 



<U9<K*)J^(o)) = o, (#Oh?(o)iaJ} = o, (13) 

where the Dirac spinor indices are omitted for simplicity. 

We want to adapt our approach and use a parametrization-free formalism, according to [l2j . To this end, we 
introduce the notations (here, T denotes different combinations of 7-matrices), 



(pa|i?(o)iv(*)bi) = ^O*) <pa$(o)r idl iK*)|pi) = rf T] {x) , 

where = denotes the Fourier transformation with the measure (z = An, z' = 0) 



(14) 



(15) 



In Eq. (fTTf , if V becomes the 7- matrix with an open Lorentz index, the functions J 7 ^ (x) and J-a F ' (x) should 
be written with additional Lorentz indices. 

Keeping the vector and axial- vector projections in the Fierz decomposition of Eq. (fT5|) (all other structures 
do not contribute in the massless quark case), the e.o.m., in terms of the functions (fl"4"|) . take the following form 



1t1~ I ?f 7+1 (s) - xP + (x) + l -e AT - a+ J^r + ^ (x) + £P+ ie?~ a+ jf T7sl (x) \ = 



7^7- j ie^ a+ T { f 7+75] (x) + (x) - xP+ ie p - a+ J^ 75 ' (x) + eP+ J*ri (x) } = 



Then, following [llj and using Eq. (fT4|) . the DVCS amplitudes (see, Eq. (ITT]) ') can be written as 



(16) 
(17) 



and 



r« = I dxix 



xP + Q 
7^7777—7777^7 



{xP + Qf 



T h+] (x) + dxtv 



XP + Q rp 

> (i_P + Q) 2>7a . 



dx tr 



xP + Q 
7^7777775777^7 75 



{xP + Qf 



jrh + ^]( x ) - dxtr 



xP + Q T 
7i/ (xP + Q) 27m7q75 



■'crossed" , 



(18) 



T« = - dx?f^(x)tr 



xP + Q T xP + Q 
7^7777— ^777 a 7-777-7777^7 



dxF [ f 7+751 (x)tr 
'crossed" . 



{xP + QY ,a (xP + Q) 2 

7m7~75 



XP + Q rp XP 

7^7777— 7777a 



(xP + Qf la (xP + Q) 2 



(19) 
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As was mentioned above, we now have to use the QCD e.o.m., written in the form of Eqs. p^|) and p7|) , for the 
amplitude T^) , and then, combining it with T$ , to collect all the similar terms in the final expression. Due 

(2) 

to the specific structure of the e.o.m., the correlators with transverse derivatives in T^J can be eliminated, and 
can be re-expressed through known correlators without derivatives (this will also be valid for the case with the 
dynamical twist contributions included). So, one gets 



T (l) + (2) = 1 

2P ■ q 



dx 



1 



+ 



1 



x — £ + ie x + !; — ie 



Tail, 



(20) 



where 



T h+] (x) 



P-M+ 7t[7 + 7b] 



J" 17 7bJ (x) + 



I — 

X 



r «-^+jr[7T7s] f x \ 



(21) 



where <5+ denotes the usual Kronecker symbol. 

We have thus derived the gauge invariant DVCS amplitude for the most general case of a target with an 
arbitrary spin which totally coincides with the results obtained in [jjj by a different approach. In the present 
paper, the corrections of the order 0(A 2 r /Q 2 ) have been neglected. We postpone the discussion of these 
corrections until a forthcoming analysis. The detailed study of these contributions has recently been presented 

m dins. 

If one now specifies spin of the hadron and, then, uses explicit parameterizations for the corresponding matrix 
elements, it will, in particular, reproduce the known cases of spin-0 and spin- 1/2 (see, [ll|, [l2j - fl7j). 

We now study DVCS off a spin-1 particle, which is of phenomenological interest in the deuteron case [(j, l23j|. 
To this end, we first specify the parametrization of the relevant matrix elements. Namely, the parameterizations 
for the vector correlators at the leading twist- 2 level are (see, 



(P2 , A 2 | [tp (Oh^(z)} tw " 2 \ Pl , A! ) £ J-]7 +1 (x) = e* 2a V$ ; J (n* , n, A T ) ei H? (x, £, t) 



(22) 



where 



4 a V^jK.n, A T )e 1/3 HY (x, £, t) = P M P^..,^, e l5 x, £, t) = 
pJ(e*-e 1 )HY(x > + (e* 2 -P)(e 1 -n) P 2 y (x, £) + (e* 2 ■ n) (e a • P) P 3 V (x, 0+ 



M 2 



(e* 2 • P)(ei • P) Hi (x, + M\e% ■ n)( ei ■ n) H£ (x,£) 



(23) 



Here, for the sake of conciseness, a new compact notation Pl\ s(e2, ei; x, £, t) has been introduced. Now, we 
are in a position to discuss the twist-3 operator matrix elements and their parameterizations. For the vector 
quark correlator we have 



(p 2 ,A 2 | $(0) 7/1 V(*)] bi,Ax> L j*rr\ x ) = £ a V^n*,n,± T )e xp GY{x,Z,t), 



(24) 



where 



4a (n*, n, A T ) ei p GY (x, i) = A'£ G^..^, e l5 x, + e^fa • P) G 6 ^ (x, + e^(e* 2 ■ P) G v 7 (x, £) + 
M 2 e\l{e x ■ n) G%(x, + M 2 e^(e| • ra) G%(x, 0- (25) 

Our next step is the parametrization of the axial-vector correlator. In contrast to the vector projection, the 
Schouten identity plays a crucial role in the determination of the Lorentz independent structures. The twist-2 
axial- vector correlator can be parameterized by 

<P2, A 2 1 [#) 7 ^W] tW " 2 |pi, Ai) £ J*r + ^ (x) = -i A { ^(n\ n, A T ) ei „ ^(x, £, t), (26) 
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where 



e 2 a 4$, fa*, n, A r )ei ^ #f (z, £, t) = ^ Pe * T ef iff (x, + ^ Va^ S - (*i • P) #2^ + 
^ 2 -^FA^(e2-P)ff 3 A (^0+^PA^e 2 ^(ei •n) J ff 4 A (x,0- (27) 

Next, let us consider the twist-3 correlators. Using the light-cone basis, we have fifteen different possible 
tensors. As in the twist- 2 case, the use of the Schouten identity reduces the number of independent tensors. 
Indeed, instead of fifteen possible structures we have nine independent tensors which parameterize the twist-3 
axial- vector correlators. Finally, the axial- vector correlator reads 

( P2 ,A 2 | [^(0)7^75^^)] tW_3 >i,Ai> = 4 TT75l (z) = -ie* 2a A%^(n*,n,A T )e ip Gf(x,£,t), (28) 

where 

4 a A%^{n\ n, A r )ei G?(x, £, t) = e^ nPeT (e* 2 ■ P) G?(x, + E/mPe . t (e x • P) G£(x, + 
M 2 £ mPe T (e* 2 ■ n) G$(x, + M 2 We . T (ei • n) Gf(x, + -^ £»A T Pe* (ex • P) G^z, + 

e M A T p e * (ei • n) Gg (x, £) + s^A T p ei (e^ ■ n) G? (x, + e^A-rne* (ei ■ P) G$(x, £) + 
M 2 s flATnei (e* 2 -n)Gi(x,C). (29) 



Inserting the explicit parameterizations ((22]) - ([28]) into the amplitude (|2T|). we derive the gauge invariant 
DVCS amplitude for the case of deuteron target: 

(\ (Ai, A 2 ) 
T« + + T$ + 7#> J + 0(A* ; M 2 ) + "crossed", (30) 

where the structure amplitudes T$y read 

r^J = Hl„ A (x; ex, 4) ^P p i>„ + P»Q» + P„Q» - g^P ■ Q) + \p»A T y - \p v A T ^j + 

GL.Ax; e u e* 2 ) ^P„A£ + 3^A^ + A^<3„ + A T V Q^ - 
Uel-P^-P) ^ + (e , p)(ei n)G A (a;)+ 

(ei • P)(e2 • n) (G$(x) - Gf(x)) ) t i^A T v - £P,A£ - A^Q U + A T V QA , (31) 



and 



and 



T$ = (ex • P)GX(x) \£P v efi + 3£P^ + e%Q v + e^Q^J + 

(ex ■ P)Gi(x) U^etl ~ ~ tfiQu + e\ T v Q A , (32) 

U 3) - (e* 2 ■ P)G 7 v (x) ^P,e^ + 3^P M eL + e^Q„ + e^Q J + 

(e* 2 ■ P)G£(x) ( $£P^ V - £P,e^ ~ ef M Q„ + ef ^ J , (33) 
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and 

T$ = £^Pn i^nPe' 2 T ef H l{ X > + Jfi e nPA^e* r ( e l ' P ) H$(x, £) + 

e nPATeI (e* 2 ■ P) -ff 3 A (x,0 + e nPATe ,T( ei ■ n) H^{x,^j . (34) 

This gauge invariant amplitude for DVCS off deuteron is our main result. For the sake of brevity, in Eqs. (|31l) 
- (|33[) . we neglected all terms which are proportional to the square of the hadron mass. The full expressions for 
all amplitudes will be presented in our forthcoming study. 



IV. SINGLE SPIN ASYMMETRY 



In the preceding section, we have obtained the gauge invariant DVCS amplitude which has a significant 
meaning for the investigation of any observables. As a phenomenologically important example, we now consider 
the single (electron) spin asymmetry (SSA) , which arises in the collision of the longitudinally polarized electron 
beams with an unpolarized hadron target. The SSA parameter is defined as 

_ d<-o - At(h 

,Al ~ <hr(-+)+d*(<-y ( ' 

The numerator of Eq. (|35|) can be expressed through the imaginary part, first, of the interference between the 
twist-2 and twist-3 helicity DVCS amplitudes and, second, of the interference between the Bethe-Heitler (BH) 
and DVCS amplitudes. 

The DVCS amplitude contributing to exclusive real photon production at Q 2 3> M 2 for the real and virtual 
photon polarizations, i and j, reads 

2 

-4dvcs = ^E L(J) A^) - ^ = C^iluhy^ > (36) 
j 

respectively. Here, the helicity amplitude is given by 

Ai,i)= e f T ^C {i) , * = ±1, 3=0, ±1- (37) 
The Bethe-Heitler amplitude reads 

4 ) H=^E A °''%)' A(i-) = V^,^ 2 )e; «el/« T -)=e(?)^, (38) 

3 

where 

A 2 = -A( 2 M 2 + A| = t . (39) 

For convenience, we introduce the following shorthand notations for Compton form factors related to the various 
GPDs: 



/ 



x — £ + ie J x — £ + ie J x x — £, + le 



and 



GYM _77v f^HV(x^) _ ^ f ^ £ G?(x,Z) 



dx 1 \ ^ =>gY, dx l \ ^ =>HY, / dx± =>g*. (41) 

a; — £ — ?e J x — — te J i i-(-i£ 

The relevant kinematics is given by 

t = -4A/V, Ay = t + 4M 2 £ 2 , M 2 = M 2 - - . (42) 
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Below, we apply the realistic forward limit \t\ <§C M 2 , such that 

Ay ~ 4M 2 £ 2 , M 2 ~ Af 2 . (43) 

We now calculate the contributions to |-4dvcs| 2 and -4-b H »4dvcs coming from the interference between dif- 
ferent virtual photon polarizations A 7 - = and A 7 * = +1 in Eqs. (j3!))) and respectively. We have the 
following expressions (here, deuteron polarizations are summed up): 

i 

[e(2&£ - i&l - uJM - m(0X Qi + Qi) + HlKW + - (i - ifW) - 

z(q£ + QXm 30Hf_e(i +oW+ a -tfW- + a - o 2 (i +0W) - 

£(2£# ~2iQX - (1 + W) + 0^ + #)((! + H)W + 

(i + 2 W + (i - O^W + W + a - 0(1 + OX) - (2£(# + #) + ^5 y ) x 

((i + 30)^r + e(i + o a ^2 - (i - 2 £^r + 0^1 - (i - 0) 2 ^D + 

(2^ + HY)(3W + CHJf - - (1 + 3£ a )ttf) + 

0^ + QX)(W + + a - OW) - OSi 4 + $X)(W - tW - a + W)} ( 44 ) 

and 



£>(°U (0)i) ] • [A (+ ^7^ +) ]* ~ { - 2G 3 (G£ + Qi + QX + # - U\ - 

i 

2£(# + SIM 2 + G 2 {-2(%X + 2£(Qi + QX))e + (2g + ixe 3 A + e 9 y )e + 
(2c - i){Qi + qX)0 + Gti-mY - m - i)(# + gX) - 

C(3C + 1)(# + ^) - 0-£i A + Qi + 6^ + 2&X + Ql - QX + U\ - U\ - 

2t{Qi Qi + Qi + QX)) + (3£ 2 + \){UX + 2£(# + QX)))} . (45) 

In the last expression, Gi,2,3 are the known electromagnetic deuteron form factors. 

As it can be seen from the parameterizations introduced above, the only surviving contributions in the forward 
limit, are related to the Compton form factors Hi, 5 and Qs^g terms. Keeping only these contributions, one can 
write 



^[L^A m ].[L^A ( 

{^ 9 y ((-i + 30< - (1 - o 2 a + oW) + ^(a + + (1 - 0(1 + 2 W) + 

nidi - e?w - (i + h 2 )W ) + - (1 + 3eiW )} + . (46) 



and 



{ - 2g 3 (qX + QX - nX)e + G 2 (m - i)GX + m + WX - ^enX) + 

Gi(-m + ^)gx - m - vox - 3< + (30 + mX)} + (47) 

All these expressions enter into the SSA. If we now calculate the imaginary part of the above-mentioned terms, 
we will obtain the numerator for experimentally accessible single spin asymmetry parameter [25| . 
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V. CONCLUSIONS 



In conclusion, we have derived the gauge invariant amplitude for the deeply virtual Compton scattering 
(DVCS) off a spin-1 hadron. As an important phenomenological application of this approach, we have considered 
the deuteron target and have presented the gauge invariant DVCS amplitude for the deuteron case. We have 
also discussed the simplest kind of asymmetries - the single spin asymmetry where the initial lepton has a 
longitudinal polarization while all other particles, the initial hadron, the final lepton and the final hadron, are 
unpolarized. 
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